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The first measurements of double-hadron production in deep- inelastic scattering within the nuclear 
medium were made with the HERMES spectrometer at HERA using a 27.6 GeV positron beam. 
By comparing data for deuterium, nitrogen, krypton and xenon nuclei, the influence of the nuclear 
medium on the ratio of double-hadron to single-hadron yields was investigated. Nuclear effects 
on the additional hadron are clearly observed, but with little or no difference among nitrogen, 
krypton or xenon, and with smaller magnitude than effects seen on previously measured single- 
hadron multiplicities. The data are compared with models based on partonic energy loss or pre- 
hadronic scattering, and with a model based on a purely absorptive treatment of the final state 
interactions. Thus, the double-hadron ratio provides an additional tool for studying modifications 
of hadronization in nuclear matter. 

PACS numbers: 13.87.Fh, 13.60.-r, 14.20.-c, 25.75. Gz 



Hadron production from a free nucleon in deep- 
inelastic scattering is generally described by fragmen- 
tation functions that contain non-perturbativc informa- 
tion about parton hadronization. These functions are 
expected to be different for nuclear targets because 
of several possible effects: energy loss of the propagating 
quarks, rescattering during the pre-hadronic formation 
process or interactions of the final-state hadrons within 
the nucleus. 

Despite recent accurate experimental data from single- 
hadron leptoproduction 0, y] and relativistic heavy-ion 



iprodu 

collisions [J, 0, @ i the underlying mechanisms in theo- 
retical models for hadronization in the nuclear medium 
differ greatly 0, S S EU El ■ In- medium modification of 
the quark fragmentation function has been described in 
terms of rescattering of gluons and quarks, and of energy 
loss due to induced gluon radiation 0, ■ Alternatively, 
colorless pre-hadron rescattering in the medium has been 
sug gest ed with additional effects due to Q 2 rescal- 
ing jllj. Older interpretations 0,0 based on hadronic 
final-state interactions require a hadron formation length 
smaller than the nuclear size, which is unlikely for struck 
quarks boosted to energies larger than a few GeV. Al- 
though models based on some of these ideas are already in 
conflict with data, clearly other types of data are needed 
to further distinguish among these interpretations. 

Double-hadron leptoproduction offers an additional 
way to study hadronization. If partonic energy loss of 
the struck quark were the only mechanism involved, it 
would be naively expected that the attenuation effect 
does not depend strongly on the number of hadrons in- 
volved, and the double-hadron to single-hadron ratio for 
a nuclear target should be only slightly dependent on 
the mass number A. On the contrary, if final hadron 
absorption were the dominant process, the requirement 
of an additional slower sub-leading hadron that is more 
strongly absorbed would suppress the two-hadron yield 
from heavier nuclei 0, ^| , so that this ratio should de- 
crease with A. 

Data from STAR [l4| on hadron pair production as 
a function of azimuthal angle showed that for a fixed 
value of the trigger hadron's transverse momentum, the 
production of opposite-side hadron pairs is completely 
suppressed for central Au+Au collisions due to the final- 



state interactions with the dense medium generated in 
such collisions. On the other hand, the same-side pairs 
exhibit jet-like correlations that arc similar to p+p and 
p+d collisions. These results were used in Ref. [l|| to 
advocate the picture that jet fragmentation occurs out- 
side the dense medium. In this model it has been shown 
that if hadron absorption or rescattering were responsi- 
ble for the observed hadron suppression, it would likely 
destroy the jet structure, and in particular the corre- 
lations between leading and sub-leading hadrons within 
the jet cone. However, the heavy-ion data cannot exclude 
hadronic absorption effects completely. 

In this paper the first measurement of double-hadron 
leptoproduction on nitrogen, krypton and xenon relative 
to deuterium is presented. All charged hadrons and w° 
mesons are considered. 

Semi-inclusive deep-inelastic scattering data are pre- 
sented in terms of the ratio 
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in which z = E^/v is the fractional hadronic energy, Eh 
is the hadron energy and v is the virtual photon energy, 
all of which are evaluated in the target rest frame. The 
values z\ and Z2 correspond to the leading (largest z) 
and sub-leading (second largest z) hadrons, respectively. 
The quantity dN Zl>0 5 is the number of events with at 
least two detected hadrons in a bin of width dz-i at z^ 
with z\ > 0.5. The quantity N Zl>0 - 5 is the number of 
events with at least one detected hadron with z\ > 0.5. 
The label A(D) indicates that the term is calculated for 
a nuclear (deuterium) target. 

The measurement was performed with the HERMES 
spectrometer using the 27.6 GeV positron beam 
stored in the HERA ring at DESY. The spectrometer 
consists of two identical halves located above and below 
the positron beam pipe. The scattered positrons and the 
produced hadrons were detected simultaneously within 
an angular acceptance of ± 170 mrad horizontally, and 
± (40 - 140) mrad vertically. 

The nuclear targets, which were internal to the 
positron storage ring, consisted of polarized or unpola- 
rized deuterium, or unpolarized high-density nitrogen, 



3 



krypton or xenon gas injected into a 40 cm long open- 
ended tubular storage cell. Target areal densities up 
to 1.4 x 10 16 nucleons/cm 2 were obtained for unpola- 
rized gas corresponding to luminosities up to 3 x 10 33 
nucleons/(cm 2 s). 

The positron trigger was formed by a coincidence 
between the signals from three scintillator hodoscope 
planes, and a lead glass calorimeter where a minimum 
energy deposit of 3.5 GeV (1.4 GeV) for unpolarizcd 
(polarized) target runs was required. The scattered 
positrons were identified using a transition-radiation de- 
tector, a scintillator pre-shower counter, and an elec- 
tromagnetic calorimeter. Scattered positrons were se- 
lected by imposing constraints on the squared four- 
momentum of the virtual photon Q 2 > 1 GeV 2 , on 
the invariant mass of the photon-nucleon system W = 
yJlMv + M 2 - Q 2 >2 GeV where M is the nucleon 
mass, and on the energy fraction of the virtual photon 
y = v/E < 0.85 where E is the beam energy. The con- 
straints on W and y are applied to exclude nucleon reso- 
nance excitations and to limit the magnitude of the ra- 
diative corrections, respectively. In addition the require- 
ment v >7 GeV was imposed to limit the kinematical 
correlations between v and z. 

Charged hadrons (i.e. 7r, K and p without distinction) 
were reconstructed for momenta above 1.4 GeV. The elec- 
tromagnetic calorimeter provided neutral pion iden- 
tification through the detection of neutral clusters origi- 
nating from two-photon decay. Each of the two clusters 
was required to have an energy E 1 > 0.8 GeV. The 7r° 
mesons were selected by requiring that the reconstructed 
invariant mass was within two standard deviations of the 
center of the n° mass peak. 

The leading hadron was selected with z\ >0.5. In this 
case, it is expected to contain the struck current quark 
with high probability. No explicit constraint was applied 
to z 2 . Both z\ and Z2 were calculated assuming that all 
hadrons have the mass of t he p ion. 

Using the code of Ref. [20j, radiative corrections to 
Rih were found to be negligible in the whole kinematic 
range. This is because there is no elastic or quasi-elastic 
tail in semi-inclusive events, and the inelastic corrections 
largely cancel in the measured ratio. 

Two methods of double-hadron event selection were 
used. Selection I contains only the combinations of 

hadron charges (lcading-sublcading) ++, , +0, 0+, 

— 0, 0—, 00. This suppresses the contributions from 
p° — > 7r + 7T - decay because the H — and — h combinations 
are missing. Moreover, in the Lund string model, the 
exclusion of the opposite-charge combinations enhances 
the rank-1 (leading) plus rank-3 (sub- leading) combina- 
tion [2lJ . The higher the particle rank, the more likely it 
is formed deep inside the nucleus, and the corresponding 
hadron absorption should be larger. Selection II contains 
all particle charge combinations. Here, the sub-leading 
hadron is mainly of rank-2 and the contribution from p° 



decay is larger. In both Selections I and II the relative 
yield from exclusive p° production in jV Zl>0 5 is small 
and was evaluated by Monte Carlo calculation to be on 
the order of 12% for the deuterium target. 
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FIG. 1: The ratio B,2h as a function of 22 for 14 N (squares), 
Kr (circles) and Xe (triangles) with z\ >0.5. Only Selection 
I was considered. The systematical uncertainty is 2% for all 
the targets and is independent of z%. In the upper panel 
the curves (solid for 14 N, dashed for Kr, dotted for Xe) are 
calculated within a BUU transport model ^(J • I n the bottom 
panel the same data are shown with calculations that assume 
only absorption for the three nuclei (same line types as in the 
upper plot) GJj]. 

Fig-fflshows the double ratio R 2 h as a function of Z2 for 
Selection I only. The kinematic variables are in the range 
{u)=2l to 16 GeV and (Q 2 )=2.1 to 2.6 GeV 2 as z 2 goes 
from 0.09 to 0.44. The averages over z 2 are (v) =17.7 
GeV and (Q 2 ) =2.4 GeV 2 . 

The ratio R 2 h is generally below unity with no signi- 
ficant difference between the three nuclei. These data 
clearly show that the nuclear effect in the double-hadron 
ratio is much smaller than for the single-hadron atte- 
nuation measured under the same kinematic conditions 
0,13- For z 2 <0.1, where R 2 h rises towards and possibly 
above 1, the slow hadrons originate largely from target 
fragmentation [uH 122^. Also for z 2 >0.4, where the two 
hadrons have similar energy, R 2 h seems to rise towards 
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FIG. 2: The ratio 7?2h as a function of 22 for 14 N (squares), 
Kr (circles) and Xe (triangles) with z\ >0.5 for Selection II. 
The systematic uncertainty is 4% (3%) for xenon and krypton 
(nitrogen) and is independent of Z2. The curves ( 14 N: solid; 
Kr: dashed; Xe: dotted) are from the parton energy loss 
model described in Refs.Hrl. ll7l . 



unity. 

Fig- n (upper panel) shows calculations based on a 
PYTHIA event generator with a fully coupled-channel 
treatment of final-state interactions by means of a BUU 
transport model 01 • I n this model, the fragmenta- 
tion function is modified by pre-hadron interactions and 
rescattering in the medium. Although the general trend 
of the data is reproduced, the model predicts an effect 
twice as large for xenon and krypton as for nitrogen above 
Z2=0.1, which is not supported by the data. 

Fig. ^ (bottom panel) shows the same data compared 
to a calculation with a purely absorptive treatment of 
the interaction of the pre-hadronic or the final hadronic 
states. The data rule out this assumption 

Fig. |3 presents Rih calculated for all hadron charge 
combinations (Selection II). Inclusion of the H — and — V 
pairs does not change the value of Rih significantly, con- 
trary to all of the naive space-time evolutionary models of 
hadronization. In order to evaluate further the possible 
influence of the exclusive and semi-inclusive p° produc- 
tion, R,2h was extracted for all hadron pairs except those 
with invariant mass near the p°. This has no noticeable 
effect on i?2ft- Therefore, the final data are presented 
over the full invariant mass range. The effect of only 
the exclusive p° production is even smaller since it con- 
tributes only 5% of the total yield. The contamination 
from exclusive production of p^ and u mesons is com- 
pletely negligible, being suppressed by more than one 
order of magnitude with respect to the p° 



The curves in Fig. [2] represent the model [lil lri| in 
which modifications of the fragmentation functions arises 
from parton energy loss. Contrary to naive expectations, 
this model predicts a significant ^4-depcndcncc, in conflict 
with the data. 

Table [I] provides a quantitative comparison between 
the data and theoretical predictions for R 2 h integrated 
over Z2- 

The total systematic uncertainty on R,2h is 4% (3%) for 
xenon and krypton (nitrogen) and is nearly independent 
of Z2- The main contribution to the systematic uncer- 
tainty comes from the decay of exclusively produced p° 
mesons. However, for the double-hadron multiplicities 
dN Zl>Q - 5 (z2)/dz2 the p° contribution has a negligible ef- 
fect. The p° contribution to jV 2l>0 - 5 was estimated in 
analogy with Ref.@, Q to be about 2% (3%) for light 
(heavy) nuclei. The only other contributing factor is the 
uncertainty in the overall efficiency of 2%. The geometric 
acceptance for semi-inclusive hadron production was ve- 
rified to be the same for both the nuclear and deuterium 
targets by studying the multiplicity ratio as a function 
of the hadron polar angle. This ratio is constant within 
experimental error. 

In conclusion, the first measurement of double-hadron 
production on deuterium, nitrogen, krypton and xenon 
is presented. The ratio of double- to singlc-hadron yields 
from nuclear targets compared to deuterium are simi- 
lar for atomic mass numbers A=1A, 84 and 131, as a 
function of the relative energy of the sub-leading hadron. 
This is at variance with the singlc-hadron attenuation 
data, which depend strongly on A. The data do not sup- 
port naive expectations for pre-hadronic and hadronic 
final-state interactions that are purely absorptive. Mo- 
dels that interpret modifications to fragmentation as be- 
ing due to pre-hadronic scattering or partonic energy loss 
are also inconsistent with the data. In fact the latter pre- 
dict an even larger A-depcndence, while the data show 
little. Like the jet correlation measurements in heavy- 
ion collisions, the double-hadron observablcs in semi- 
inclusive deep inelastic scattering provide new informa- 
tion for differentiating between models of hadronization 
in nuclei that are indistinguishable in singlc-hadron mea- 
surements. 
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TABLE I: Averaged values of R 2 h for Z\ >0.5. 





22 <0.5 


0.K 22 <0.5 


0.K 22 <0.5 
Model [10] 


0.K 22 <0.5 
Model Abs. [10] 


0.K 2 2 <0.5 
Model [16. 17] 


i4 N/D 
Selection I 


0.946 ± 0.017 ± 0.019 


0.941 ± 0.018 ± 0.019 


0.931 


0.907 




i4 N/D 
Selection II 


0.975 ± 0.009 ± 0.029 


0.972 ± 0.010 ± 0.029 






0.965 


Kr/D 
S6lection I 


n qoq 4- n rn ^ -i- n mo 
u.yzy ± u.uio zt u.uiy 


u.yi* it u.uio =t u.uio 


U.ooO 


n 7Q^ 
u. / yo 




Kr/D 
Selection II 


0.902 ± 0.008 ± 0.036 


0.892 ± 0.008 ± 0.036 






0.879 


Xe/D 
Selection I 


0.936 ± 0.023 ± 0.019 


0.915 ± 0.024 ± 0.018 


0.815 


0.773 




Xe/D 
Selection II 


0.936 ± 0.012 ± 0.037 


0.925 ± 0.013 ± 0.037 






0.800 
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